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I. INTRODUCTION 
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We calculate the frequency dependence of gravitational wave background arising at second order 
of cosmological perturbation theory due to mixing of tensor and scalar modes. The calculation of 
the induced gravitational background is performed for two special cases: for the power spectrum of 
scalar perturbations which has a peak at some scale and for the scalar spectrum predicted by the 
inflationary model with the running mass potential. We show that the amplitudes of the induced 
\ gravitational background, in the frequency region ~ 10~^ — 10'' Hz, are effectively constrained by 

. results of studies of primordial black hole production in early universe. We argue that though 

' today's LIGO bound on Vn(k) is weaker than the PBH one. Advanced LIGO will be able to set a 

^Nj \ stronger bound, and in future the ground-based interferometers of LIGO type will be suitable for 

^ . obtaining constraints on PBH number density in the mass range ~ 10^^ — 10^^ g. 
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u: According to the one of the main predictions of General Relativity, the present Universe is filled with a difi^use 
(-H ' gravitational background produced by sources of astrophysical and cosmological origin (see reviews [H, Q). In 
O ^, particular, relic stochastic background necessarily arises in all standard and superstring-motivated ("pre-big bang" 
' (PBB) i4:]) inflationary scenarios, as a result of the process of the amplification of vacuum fluctuations (such an 
2 , amplification of quantum fluctuations of the gravitational field modes had been discovered in [Bj and had been used 
in a model with de Sitter phase of expansion in [^). The standard inflationary scenario is characterized by a flat or 
^ slightly decreasing spectrum 0,B&lld|j which is constrained at the scale of the present Hubble radius (at / ~ 10"^^ 
Hz) by the large scale isotropy of the CMB radiation (Til . [T^ . On the contrary, PBB models predict a growing 
spectrum [l^, [3| constrained at high frequencies by the nucleosynthesis bound [3, [lB| ■ 
^ , In the period near the end of inflation and after it the processes are possible which could result in additional 
' gravitational background: in particular, first order phase transitions llTI . |l8| (such transitions can occur even before 
\Q , the end of inflation, like, e.g., in some two-field inflationary models [13|) and preheating (see, e.g., [lOiHH). These 
processes lead to gravitational wave radiation due to local strong inhomogeneities generated in the cosmological fluid 
, during the transition from inflationary expansion to radiation era and subsequent phase transitions. The corresponding 
00 ' spectra of the background radiation have a peak at some frequency (the peak's position depends on the inflationary 

' energy scale and on the temperature of the phase transition), see, e.g., (22l. f23|. 
0^ . The rather substantial gravitational wave background (GWB) is possible in cosmic string scenarios (see, e.g., the 
' recent work (24|). in brane world models (see, e.g., p^). in scenarios with intense production of p rimordial black holes 
J> (PBHs) in early universe (in the latter case GWB appears as a result of PBH evaporations [la[13|)- 

It was realized during last decade that there is still another background of stochastic GWs of cosmological ori- 
gin. Namely, GWB arises as a result of non-linear evolution (i.e., of gravitational instability) of curvature (density) 



■ fluctuations. The nonlinear effects (the mixing of different modes) appear already in second-order of cosmologi- 
cal perturbation theory [2^, [2^[33|- In particular, a second order contribution to the tensor mode, h^i^j' , depends 
quadratically on the first order scalar metric perturbation, i.e., the observed scalar spectrum sources the generation 
of secondary tensor modes. By other words, the stochastic spectrum of second order GWs is induced by the first 
order scalar perturbations. Calculations of Q.aw at second order and discussions on perspectives of measurements of 
the second order GWs are contained in works [3l|, [13, [H, H^l ■ 

The generation of GWs from primordial density perturbations on very small scales which are not directly studied 
by astronomical measurements could be used for constraining overdensities on these scales, in a close analogy with 
the case of primordial black holes [H, [s^. Large curvature perturbations leading to "features" (e.g., peaks or spikes) 
in the primordial power spectrum and to possible PBH production, can arise in multiple field scenarios at the end 
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of inflation (during the preheating era) or between two consecutive stages of inflation, as a result of parametric 
resonance or tachyonic instabihty. Such features can, in principle, exist even in single-field inflationary models (see, 
e.g., [3^). These peculiarities of the primordial scalar spectrum lead to local enhancements in the induced spectrum 
of gravitational perturbations [3^ |34| . Another example of an inflationary model predicting large amplitudes of the 
density perturbations at small scales is the inflationary model with the running mass potential. In a case of the 
positive running, the scalar spectrum at large values of k (near kend) can steeply rise with k (which is also a kind of 
the "feature"). 

The aim of the present paper is two-fold: we calculate a spectrum of induced GWs for two cases, which are not 
studied in previous works [3^. l33l. [33 | : for a scalar power spectrum with a peak of non-zero width and for a scalar power 
spectrum with the running of the spectral index (the latter case is studied with using the particular model, namely, 
a model with the running mass potential [37l . [ssjV The second aim of the paper is a constraining of the induced 
GWB using the results of PBH searches. We consider, in the present paper, the tensor spectrum in a rather narrow 
interval of wave numbers, corresponding to modes leaving horizon at the time near the end of inflation. These modes 
enter the Hubble scale during the radiation-dominated (RD) era. Overdensities lead to production of PBHs with 
small masses (~ 10"'^^ — 10"'^^ g) which have enough of time for evaporation during the life of the universe. Products of 
evaporation of these PBHs contribute to extragalactic diffuse photon and neutrino backgrounds (which are measured 
experimentally). This allows to obtain constraints on the primordial power spectrum amplitudes. An independent 
constraint on the scalar spectrum can, in principle, be obtained by a direct detection of induced GWs [s^jUHl- In this 
paper we compare the abilities of two methods of such a constraining. 

The plan of the paper is as follows. In Sec. [H] calculation of the induced GW background is performed for 
two different cases: for the delta-function power spectrum of the primordial scalar perturbations and for the power 
spectrum with a peak of finite width. The comparison of PBH constraints for the induced GWB with possible 
constraints from the future experiments, such as Advanced LIGO, is given. In Sec. Illll the analogous calculation and 
comparison are done for the scalar spectrum predicted by the running mass model. The last section contains our 
summary and conclusions. 

II. GW BACKGROUND CALCULATION 

A. Connection between frequency and horizon mass 

For a wave with comoving wave number k and wavelength A = 27r/fc, propagating at the speed of light c, the 
corresponding frequency is / = c/A, or 

/ = £^ = 1.54 X 10-1^ f—^^ Hz. (1) 
2n VMpc-V 

From the constancy of the entropy in the comoving volume, we have the relation between the scale factor a, temper- 
ature T and the effective number of degrees of freedom g„: 

a^g^^'^'T-K (2) 

From the Friedmann equation {H^ ^ p), we have 

Hr^a-'g-'^', (3) 
and the horizon mass corresponding to the scale factor a evolves during the radiation-dominated (RD) epoch as 

Mh^iH-'fp^a^g'J''. (4) 

From ^ and ([?]), the wave number of the mode entering horizon at the moment of time t (at this time, k ~ aH) is 
related to the horizon mass at the same moment of time by 

^ ^ (^) (£;) " ^ 10^^(M,[g])-V2 Mpc-\ (5) 
where in the last equality we have adopted that g^^^q ~ 3, w 100, 

Meg = 1.3 X lO^^g • {n^h^^ « 8 X lO^^g, (6) 
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keq = aeqHeq = V2Honmi^j^^^^ « 0.0095 Mpc'^ (7) 

The frequency of the wave corresponding to the wave number k can be related to the horizon mass by the relation 
following from ([T]) and ([S]), 

Q X 1 0^^p- 

/ « 3 X 10« Hz X (M,[g]ri/2. (/[Hz])2 - 

For scalar-induced GWs, the single mode in scalar spectrum does not correspond to the only one mode in Vh- For 
example, for the (5-function-like spectrum V-jzik) ^ S{k — fco), the GW spectrum is continuous and stretches from 
to 2ko [SJ]- However, the order of magnitude of wave numbers of induced GWs, as we will see, is the same as of 
scalar perturbations, so ([8]) gives an estimate of GW frequency that will be generated from perturbations entering 
horizon at its mass scale Mh- Furthermore, if PBHs form from a scalar spectrum of perturbations at a horizon mass 
scale Mh, the typical PBH mass will be of order of Mh (see, e.g., [i^l), so ([8]) relates the typical PBH mass with the 
characteristic frequency of second-order GWs produced. 



B. General formulas 



The components of the spatially flat FRW metric are following: 

500 = -a2(l + 2$), 
{diUj + uji) , 



go 



9ij 



(1 - 2^-) (5y + D,jh + d,hj + djhi + -K^ 



(9) 
(10) 

(11) 



For a derivation of the equation of motion for second order tensor perturbations we used the generalized longitudinal 
gauge, which is defined by the relations 



= /l(l) = h. 



(1) 



0. 



(12) 



Besides, we neglect, following [ssj, the first order vector and tensor perturbations. 



(13) 



i(2) 



With these approximations the components of the second order Einstein tensor, , depend only on first-order 

perturbations <i>(^\ ^!^^^ and second-order perturbations ^'^^^ uj^'^\ hP'K h!"^^ and /i^^"*. The second-order 

energy-momentum tensor, with ignoring anisotropic stress contributions (in this case $ = ^'), is given by 



t: 



(2) 1 



(1) 



(14) 



Here, u^s is the vector (solenoidal) part of the first-order velocity perturbation. 



3^2(1 + u;) 



(15) 



where Ti = aH and w = p'^'/p^"^ is the equation of state parameter. Using the Einstein equations, with keeping in 
them only transverse-traceless contributions, 



G: 



,j(2) 



TT 



SttG 



J (2) 



TT 



(16) 



leads to the evolution equation for the second-order tensor perturbation hij , which contains, in the source term, only 
first-order scalar perturbations (all the second-order perturbations except of h\-' are eliminated): 



h'lj + 2Hh'ij - W^h.j = -4 



TT 



(17) 
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Here, we have omitted the indexes: "^^^^ 5', h[^'' hij. The ehmination of all second-order perturbations except 
of h)^ is due to using of transverse traceless projection in Einstein equations (see Eq. (|16p ). It had been shown in 

[33 , Isst (and can be straightforwardly checked) that all the second-order terms which enter the Einstein tensor (except 
(2) 

of hlj) are canceled after such a projection. 

According to [ssj , the power spectrum of induced GWs is given by the expression 

00 1 

Vh{k,T) = j dk j dfiV^i\k-k\)V^{k)^{k,k,^i,T), (18) 
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where 



J^{k,k,fi,T) = „, , ~- / dTi a(fi)5fc(T,fi)/(k,k,fi) x 

To 

r 

X y"df2a(T2)5fc(r,f2) [/(k,k,f2) + /(k,k-k,f2)] (19) 



To 



and 



/(k, k, t) = 12*(fcT)*(|k - k|r) + 8T^ikT)^'(\k - k|r) + 4r2*'(fcr)^''(|k - kjr). (20) 



In Eqs. UHl Uni HO]) the following notations are used. V^{k) is the power spectrum of the Bardeen potential, defined 
at some moment of time r — near the beginning of the RD stage (by definition, it is the primordial spectrum). 



^I^k is the Fourier component of '5, 



(*k*k') = ^'53(k-fk')7'*(fc), (21) 



*W = / d'^^e^-'^ (22) 



/i = k • k/{kk) is the cosine of the angle between the vectors k and k. The power spectrum of GWs is defined by the 
standard expression 

{h^{T)h^'{T)) - ^'^6Hk + k')Vh{k,T), (23) 
where h^^ir) is the Fourier component of the tensor metric perturbation, 

h^A^, r)^J (^^e^"^-" [hu{T)e^,{k) + h^{T)e,,{k)] , (24) 

Cij (k) and (k) are two polarization tensors corresponding to the wave number k. 

It had been shown in works of previous authors [3^ . [33| that the right-hand side of this equation is expressed 
through the correlator containing the product of four \l/fc functions, i.e., Vh is expressed through the product of two 
V\si spectra, see Eq. (18). So, in this case, like in the case of Eq. ipT]) . the average is also over fluctuations of the 
gravitational potential. 

The evolution equation for the GW amplitude is 

K + + k^hk = S{k, r), (25) 

where the source term is 

^(k, r) = 1 d'k P(l - ^,^) /(k, k, T)^k_i,^i,. (26) 
The function / in Eq. (|20p contains transfer functions ^'(fcr), which are defined by 

*(fcr) = -^, (27) 
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where = ^'a:(Tj') is the initial (primordial) value of the potential. During RD epoch, the solution for the Bardeen 
potential, havin g th e initial condition ^'^(Tj) = 0, where is the moment of the end of inflation which is close to r^' 
(but n < tI), is |3i,|43] 

^kir^TZk)^ — :t-\{x — Xi) cosix — Xi) — {1 + xxi) su\{x — Xi)\ , x — kr/V?,. (28) 

X'^ 

Here, the variable TZ is the curvature perturbation on the conioving hypersurfaces (see, e.g., (39l|). 

The value of the potential at is chosen to be zero because is typically very small during inflation [SO^ and 
it is a continuous function during the transition from inflationary to RD stage (we assume, for simplicity, that the 
reheating is instant). For matter-dominated (MD) epoch, ^'fc(T) = const on all scales. 

We have chosen r^' using the condition lg(Tj'/Tj) = 0.05, and due to this the Bardeen potential '^^ at t'^ is much smaller 
than its asymptotic value = —(2/3)7?.^ that is reached in the super-horizon regime {k <C aH) for k <C k^nd = 
(this result can be obtained from (|28p expanding the sine and cosine functions). If we are interested only in such wave 
numbers (fc <^ kend), it is more convenient to deflne V^i in terms of this asymptotic super- horizon value. To distinguish 
asymptotic value of the \l/-spectrum from the value at the moment r/, we will denote it as P*. For k <C kend, the 
relation between the two is very simple: 

and 'J'fc — ^'fe(Tj', 7?.fe). For our choice of t[, k w 11.7. 

The function gk (t, f) in Eq. is the Green function of the Eq. (pS)) which depends on the cosmological epoch. 
For RD Universe, 

9k{T, t) = ^ sin[fc(T - •?)] , T < Teq, (30) 

and for MD case, 

XX 

9kiT,f) = -—[ji{x)yi{x) - ji{x)yi{x)], x = kr, r > Teq. (31) 

C. Delta function input power spectrum 

The integral (|18[) is much simplified if we assume an idealized power spectrum with all power contained in one mode 
with some wave number /sq: 

P*(fc) = Pq5 {\n = PokoS{k ~ ko). (32) 

Such a case has already been studied in [H, [13]. The formula for the GW power spectrum for such an input is 
obtained from ^5)) : 

P/.(fc) = ^^(^A:,fco,M=^,rj. (33) 

From /i — cos 6 < 1, it follows that k < 2ko, and GWs in this case are generated in the frequency interval from to 
2fco. 

The corresponding asyptotic value of the coefficient in the delta-spectrum ((32)) is equal to Pq — kPq. 

D. GW energy density calculation 

The energy density of GWs per logarithmic interval of k in units of the critical density is given by 

I / k 



12^^1(7)17(7)' ^''^ 



Ig a 



FIG. 1: The dependence of Vh{k) on the scale factor a for several wave numbers. As an input, we used here a delta-function 
power spectrum for Vm, with Pq = 10""^ for each case. For curves from top to bottom, k = ko — 6keq, 20keq, 3 x lO^keq. 

The power spectrum of GWs Vhik, r) is obtained from the formula (llSp . However, for very large wave numbers k which 
we are interested in, the direct use of (|18p will require numerical integration for functions having a huge number of 
oscillations (e.g., for k ^ lO^'^^Mpc^^ this is about ^ fcro ^ 10^° oscillations). This is hard to do numerically. 
Fortunately, we do not have to do integration until the present day. It is enough to calculate fJgw for the moment of 
time Tcaic ^ at which the mode is well inside the horizon, and is freely propagating. We can then easily relate 
energy densities of GWs at different times with simple calculation, using the fact that hk ^ far inside the horizon. 
Really, assuming zero source term in (j25p and changing the variable to Vk — ahk, we obtain the equation 



where w = p/p. In sub-horizon regime, when k 3> aH, the solution of this equation is Vk ~ cos(fcT -I- ip), and so, 
ignoring oscillations, hk ^ . 

However, it had been noticed in 33] that the propagation of second-order GWs in sub-horizon regime cannot be 
always regarded as free because the source term in the equation (|25p cannot be neglected in all cases. We illustrate 
this point in Fig. [Jl which shows the dependence oi'Ph{a) for several values of k, calculated numerically using Eqs. 
(fT51 [501 ED). The input power spectrum was taken to be of a delta- function form for this example. It is seen from 
this figure that for rather large values of fc (fc > fee ~ lOOfceij ~ lMpc~^) this effect can be neglected, and a simple 
relation Vh{a) ~ can be used while for smaller k the effect is important. In the present work we are interested in 
GWs with much larger wave numbers than fee, so we can safely use this relation in our calculations. 

During the RD epoch, aH ~ a~^g7^^^, so i7gw(fc,T) ~ {k/aH)^T'h ~ gV^, and we can write the relation for 
moments of calculation Tcaic and matter-radiation equality T^q 



After the moment r^q, f^cw is proportional to energy density fraction of the radiation, which equals 0.5 at r^q. So, 




(35) 




(36) 



0° 



'gw(^) — 251 7? X Q^^lk) 



(37) 



and, finally. 




(38) 



This formula gives the correct energy density, accurate to the oscillations in it. The exact shape of the function will, 
actually, depend on the choice of Tcaic 3> k~^, and the larger Tcaic we take, the more frequent are the oscillations, but 
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14.5 15.0 15.5 16.0 16.5 17.0 

lgk[Mpc"^] 

FIG. 2: GW spectrum from a delta-function peak in Vi {Po = 1-2 x 10~^, feo = 8 x 10^*^ Mpc~^). Thin solid line - calculation 
for Tcaic ~ lO^k^^ , thin dashed line - for Tcaic = (50 + Nrud)k~^ , thick line is the envelope. 

the envelope which we are interested in does not change. In practice, Tcaic can be either fixed or dependent on k, e.g., 
for the last case, 

Tcaic = iVsub • TVsub-lOO. (39) 

It proves to be more convenient to use the "randomized" value of Nguh, i-e-, 

Tcaic = (^sub+A^rnd)-fc~\ (40) 

where A'sub is constant and A^rnd is a random number in the interval [0, 27r] calculated independently for every k. In 
this case the result of the calculation is a stochastically oscillating function whose envelope always can be easily found, 
and it is the envelope that we are interested in. This argumentation is illustrated in Fig. [2] where flaw is calculated 
for the delta- function power spectrum using two different choices of Tcaic (constant and dependent on k). It is seen 
that the envelope is the same for both cases. 

It should be noticed here that for an input spectrum of an ideal (5- function form, the growth of Vhik) for k — 2ko /VS 
proceeds logarithmically even when fcT ^ 1 |32| . It happens because of the "resonance" between the functions g and 
/ in integral in the RD epoch. Indeed, omitting the constant phase shifts, 

fffc(T,f) ~ sin(fcT); (41) 



1 . I kr \ . / |k — k|T \ ^ ■ 2 f ^qt 



_Ui.^ ._„„^(^^|, ,42, 



and amplification during integration is possible if the condition 

holds, i.e., for the case k = 2fco/\/3- The width of the resonant peak around this value of k is proportional to {kT)~^ 
and its height ^ \n{kT) [s^l, so the power contained in the peak is small and hardly can be detected. For a realistic 
spectrum of a finite width, the "resonant" effect still exists, but the amplification continues only until t ~ 1/ Ak, with 
Afc being the characteristic width of the spectrum. 
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FIG. 3: Calculation of Qowik) at the present epoch for finite width curvature perturbation power spectra of the form (|44| 
(from bottom to top, E = 0.1,0.3,0.8 and = 0.01; upper curve is for scale-invariant input spectrum with Vnik) — 0.01) 
We assumed that (?,(fco) ~ 100. 



E. Power spectrum with maximum 



It is convenient to use some kind of paramctrization to model the realistic peaked power spectrum of finite width. 
We use the distribution of the form 



IgPK(fc) =B + (lgP^-S)exp 



(Igk/ko) 



(44) 



Here, B « —8.6, characterizes the height of the peak, ko is the position of the maximum and S is the peak's 
width. Parameters of such a distribution have been constrained previously 40] from non-observation of PBHs and 
products of their Hawking evaporation. 

In Fig. [3] we show the result of ficiv calculation for the finite- width distribution of the form (|44|) . It is seen 
that for a narrow peak, the distribution looks hke the one produced by a ^-function power spectrum. The shape is 
smoothing with the growth of S and it is scale- invariant for the scale- invariant input. The value of ilcw in this case 
is proportional to [V-jiY and can be estimated as 



VtGw{k > fcc,ro) ^ 0.002 



9*eq 
9* 



1/3 



■ n 



(45) 



It follows from this formula that for obtaining large second-order GW background we need relatively large value of 
Vu- For instance, for VLqw to be of order of 10^'', Vu should be of order of 10~^. In our paper we are considering 
examples in which the values of Vn are just so large on small scales. The only constraint on V-r for these scales comes 
from PBH studies, and the values can be, in principle, much larger than V-ji ^ 10^^ obtained from experiments at 
large scales. 

Note that the shape of the curves shown in Fig. [3] does only depend on the relation fc/feg if other parameters are 
fixed and k, fco ^ kend- It can be seen from (flSl [T9|) : if we make a shift of k (and k) by a constant factor of a (so, 



k — > ak, k — > afc, and V^{k) V\s/{ak)), and the corresponding shift of calculation time Tcaic — > a 
kr^aic remains unchanged), then functions entering the Eq. pop change in the following proportion: 



g~ isin(fcr)~ i; f ^ ^^kr) ^ I- T^k^drfg^f 
K a 



1 

a 



Tca/c (so that 



(46) 



(the change of the integration limits over r does not affect the result as soon as tq <C fc ^). Using (|46)l and p8|) . we 
obtain that 



Vhik^Tcalc) 



dk / dfiV^i,VM,J' ~ a° 



(47) 
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FIG. 4: Calculation of Qcw{k) at the present epoch for the case fo = 100 Hz, E = 
— 0.016 (lower curve). Such parameters are maximal allowed from PBH constraints, 
on Q.GW obtained in the LIGO experiment and bound range expected in the future. 



3, Vji = 0.032 (upper curve) and 
Also shown are experimental limits 



SO Vh, and, correspondingly, Qqw do not change at the calculation time, and, in fact, will have the same values after 
the shift of k because, as we have seen, the amplitude of the real GW spectrum does not depend on Tcaic (as soon as 

Tcalc > k-^)- 



F. Comparison with current experimental data 

The ground-based interferometer LIGO during its fifth science run (S5) have obtained the limit [4l| 

new < 6.9 X 10"^ (48) 

This limit applies to a scale-invariant GW spectrum in the frequency range 41.5 — 169.25 Hz. The previous limit (S4 
result) was about an order of magnitude higher [4^ . 

naw < 6.5 X 10-^ (49) 

for the frequency range 51 — 150 Hz. The target sensitivity of the planned Advanced LIGO experiment is Haw ~ 

10-® - 10-9 m. 

The corresponding horizon mass for the central frequency of LIGO sensitivity range / ~ 100 Hz calculated from 
([5]) is about 10^^ g. From the other side, it is known that the intensity of PBH production in early universe from 
scalar power distribution of the form (|44p in this mass range can be constrained from studies of photon and neutrino 
extragalactic diffuse backgrounds [i^. These constraints lead to corresponding limits of V^. Particularly, for peak 
with horizon mass corresponding to its maximum M° = 10^'^ g and width E = 3, the constraints obtained in j^^l 
are V^f^ = 0.016 (in case of standard collapse model, for which critical density contrast leading to PBH formation 
is Sc = 1/3) or — 0.032 (for the critical collapse with adopted threshold value Sc = 0.45). The most optimistic 
allowed signal of second-order GWs in the LIGO range can thus be estimated from (|45|) : 

Qqw «^ 0.002 X (3/100)^/^ X 0.032^ « 6 x 10"^ (50) 

which is an order of magnitude smaller than the current bound pS)) . but is reachable for Advanced LIGO. The results 
of the full calculation of Oci^-distribution expected in this case are shown in Fig. 21 It is seen that Advanced LIGO 
will be capable to reach sensitivity needed to improve limits on Vn and PBH abundance in this range of scales. 
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Ig k [Mpc '] 

FIG. 5: Power spectrum Vnik), calculated for the running mass model, with Trh = 10^" GeV and no = 0.96; nj, = 4.5 x 10"^ 
for the upper curve and no = 4.0 x 10""^ for the lower curve. 

III. RUNNING MASS MODEL 

The runninff mass inflation model was proposed in [stI [s^ and further studied in many papers including j43l . 0, |45| . 
m, 53 and [33 . The model predicts a rather strong scale dependence of the spectral index, possibly allowing large 
values of V-jzik) at small scales, which can even lead to significant PBH production and helps to constrain possible 
model parameters [S, H Ei, [sq] . 

The potential of the running mass model takes into account quantum corrections in the context of the softly broken 
global supersymmetry and is given by the formula 



-TO^(ln 



(51) 



Its shape is determined by parameters c and s, which are defined by the relations 



'M2 



dm 



dim 



0=00 



s = c In — 
4>o 



(52) 
(53) 



Here, 0, is the inflaton field value corresponding to the maximum of the potential, 0o is the value at the epoch 
of horizon exit for the pivot scale fco ~ 0.002/i Mpc^^. The same parameters determine the behavior of the power 
spectrum of density perturbations. The connection of s, c with the observable quantities no and Uq (the spectral 
index and its running, respectively) is given by 



no - 1 « 2(s - c) , no « 
According to analysis of [s^, the possible choice is 

c« 0.06 ; s « 0.04, 



2sc . 



(54) 



(55) 



corresponding to no ~ 0.96 and u'q « 0.005. In the case of positive s and c, the inflaton field decreases during inflation. 
It means that the slow-roll parameter e also decreases: 



V^Vo 



(56) 



Correspondingly, the behaviors of the power spectra for density perturbations and gravitational waves are strongly 
different, 



Vo (infl) 



Vo 



(57) 
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FIG. 6: Calculation of Qcwik) for the power spectrum of curvature perturbations expected in the running mass model (see 
Fig. [5]). We used the following set of parameters for the calculation: Trh = 10^'^ GeV, no = 0.96 (both curves); uq = 0.0045 
for the upper curve and Uq = 0.004 for the lower one. Also shown are the experimental limits on ^Igw obtained in the LIGO 
experiment and bound range expected in the future. 



As a result, the power spectrum of induced GWs may be quite substantial at small scales in spite of the fact that 
T^Ij'"^^ is negligibly small. The parameter s connects the field value (/)o with the Hubble parameter during inflation 
and with the normalization of the CMB power spectrum: 

cl^os - 5^ . (58) 

From a theoretical point of view [i^l, Hj can lie in the wide range of values, depending on the mechanism for 
supersymmetry breaking, from Hj ~ 10"* GeV for "anomaly-mediation" case to Hj ~ lO^'^ GeV for "gauge-mediation" 
(we suppose that the inflationary potential is of order of Mf^^i where Minfi is the scale of supersymmetry breaking 
during inflation which is approximately equal to the scale of supersymmetry breaking in the vacuum). Strictly 
speaking, the form and amplitude of the power spectrum depend on the value of Hi, according to Eq. (j58p . But 
we will assume, for simplicity, that the power spectrum behavior is determined mainly by parameters s and c, 
independently on Hi, while the value of Hi determines only the temperature of reheating, Tmii and a value of the 
comoving scale crossing the horizon at the end of inflation, k^nd- If the reheating is instant, these values are given by 



^end 



Mpc"' . (60) 



The tensor power spectrum generated during inflation and the amplitude of the inflationary GWB are determined 
by the energy scale of inflation 



(61) 

k=aH 



and, because H « Hi, we have P^'"''^ < 10-30 even for the largest possible value of Hi. This is too small to be 
detected even for the most sensitive experiments proposed. 

The maximum values of frequencies of the GW spectrum depend on the moment of time of phase transition to RD 
stage from the preceding stage (see, e.g., Q). For more general discussion about the form of the relic GW spectrum 
see the recent review '5ll|. In particular, in inflationary models, the maximum value of the frequency depends on 
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Igf [Hz] 
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-6 
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-14 
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Ig k [Mpc '] 

FIG. 7: Calculation of i}Gw{k) for the running mass model. For this Figure, we used the following set of parameters in the 
calculation: Trh = 10* GeV, no — 0.96 (both curves); uq = 0.0054 for the upper curve and ng = 0.005 for the lower one. We 
also show the expected sensitivity curve for the proposed BBO experiment [5^]. 

the reheating temperature, see, e.g., Eq. (j60p . Our model with rumiing mass potential predicts large values of the 
curvature perturbation spectrum just near kgndi i-e., the position of the maximum of second-order GW background 
also depends on the reheating temperature. 

The numerically calculated power spectrum V'jz{k) in the running mass model is shown in Fig. [5] for two sets of 
parameters. Parameters uq, n^of the spectra are chosen with taking into account the PBH constraints obtained in 
the previous work of authors [40|. It is seen from the figure, that the maximum values of the power spectrum (at large 
k) are very sensitive to the value of the spectral index running. The upper curve in Fig. [5] represents the spectrum 
with largest possible for other parameters being fixed. Larger values of spectral index running will cause too large 
curvature perturbation spectrum values - in this case, too many PBHs will be produced with mass ^ 3 x 10^^ g. 

The spectra of second-order GWs corresponding to the scalar spectra of the running mass model are shown in Figs, 
m [T] As in the previous case of a peaked-power spectrum (Fig. |4]) , it is seen from these figures that scalar perturbation 
power spectrum generated in the running mass model can be a source of significant amount of GWs, detectable in 
experiments studying the frequency region ^ 10~^ — 10^ Hz. 

The results shown in Figs. [6l [7] correspond to different values of kend and Tjih (i-c, to different energy scales of 
inflation). It was shown in [4^ that PBH constraints depend on a proposed value of Tj^h'- the maximum running 
allowed is nf, = 0.0045 for Trh = 10^° GeV and = 0.0054 for Trh = 10* GeV (for hq = 0.96 in both cases). The 
difference in PBH constraints originates from the connection between comoving scale and horizon mass (see Eq. [SJ: 
k ~ kend ~ 6 X 10^^ Mpc~^ corresponds to Mbh Mh ~ 3 x 10^^ g, while k ~ k^nd ~ 6 x 10^'' Mpc~^ corresponds 
to Mbh ^ Mh « 3 x 10"'^"'^ g. PBH constraints in the region of Mh ~ 10^^ g follow from inspection of extragalactic 
photon background data while PBH constraints near Mh ~ 10^^ g follow from available limits on the intensity of 
extragalactic neutrino background |40| . 

IV. CONCLUSIONS 

The main new results of the paper are shown in Figs. [31 |3|and[Sl [71 

1. The dependence of a form of the ^lcw{k)-cva\^e on the width of the peak of the scalar power spectrum is 
carefully studied (Fig. [3]). It is shown that this curve has a characteristic "double-peak" form; the height of the larger 
maximum weakly depends on the width of the proposed Vnik) distribution (see Eq. (jll])). It is shown also that the 
0(3 vK(fc) -function depends, in fact, on the ratio fc/fco where k^ is the position of the maximum in V-jzik)- 

2. The distribution flowik) for the induced GWB is calculated for two different cases: for the scalar spectrum 
V-iz{k) with a peak (Eq. (j44p ) and for the scalar power spectrum predicted by the running mass model. It is seen 
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from the resulting figures that the behavior of {k)-cnTYes near the maximum in both cases is rather similar. 

3. It is shown that maximum values of r2GH'(^)"distributions for the induced GWB can be constrained by PBH 
searches in the wide region of the comoving scales, 10^'* ^ k ^ 10^* Mpc~^ (the upper curves in Figs HI El [7] are 
maximum ones in this sense). It is shown also that, up to now, the PBH constraints are more strong than the 
available limits from LIGO (see Fig. [4] and [6]). In particular, it follows from Fig. |4]that the maximum value of flow 
at frequency / ^ 100 Hz is equal to 6 x 10^ (supposing that the scalar power spectrum has the peak in this region). 

Finally, it follows from Fig. 3] that, in not very distant future, when the limits for Qgw from LIGO experiment will 
reach the level ~ 10~^ and lower, their data will give stronger constraints on amplitudes of scalar power spectra then 
today's constraints following from PBH studies. It means, that the future experiments with LIGO detector, limiting 
the power spectrum of primordial scalar perturbations, will give the new constraints for PBH production in the early 
universe. The interval of frequencies of ground-based detectors is about ~ 10 Hz to ^ few kHz, and the corresponding 
interval of PBH masses is from ~ 10^^ to ~ 10^^ g. 

In this paper we have made an accent on second-order GWs with frequencies ~ mHz - kHz, which will be probed 
by space-based and ground-based laser interferometer experiments. The same effects (particularly, rather large values 
of induced flow) can appear in other frequency regions, for which Viz is not probed by observations (/ ^ 10~^^ Hz), 
e.g., for / > 10^ Hz, where high-frequency GW detectors operate, or at / ~ 10~* Hz, where ilcw is constrained by 
pulsar timing data (in fact, these data already allow to put limits on PBH abundance [s^] for PBHs with rather large 
masses). In the band probed by CMB measurements, the effects from second-order GWs are carefully studied in |33j 
and shown to be small. 
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